Fixed Point Properties of the Ising Ferromagnet on the Hanoi Networks 
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The Ising model with ferromagnetic couplings on the Hanoi networks is analyzed with an exact 
renormalization group. In particular, the fixed-points are determined and the renormalization- 
group flow for certain initial conditions is analyzed. Hanoi networks combine a one-dimensional 
lattice structure with a hierarchy of long-range bonds to create a mix of geometric and small-world 
properties. Generically, those small-world bonds result in non-universal behavior, i.e. fixed points 
and scaling exponents that depend on temperature and the initial choice of coupling strengths. It 
is shown that a diversity of different behaviors can be observed with seemingly small changes in 
the structure of the networks. Defining interpolating families of such networks, we find tunable 
transitions between regimes with power-law and certain essential singularities in the critical scaling 
of the correlation length. These are similar to the so-called inverted Berezinskii-Kosterlitz-Thouless 
transition previously observed only in scale-free or dense networks. 



INTRODUCTION 



The study of equilibrium statistical models, in partic- 
ular the Ising model, on complex networks have led to 
a number of novel phenomena, many of which still defy 
proper classification 0, The reasons for that orig- 
inate with the diversity of network structures conceiv- 
able, where the most basic features, such as their degree 
distribution or average range of connections, are insuf- 
ficient to characterize their asymptotic properties. In- 
stead, even within a given ensemble of networks, scaling 
behavior may depend strongly on the details of the net- 
work structure or the strength of the interactions. De- 
spite of that high degree of non-universality, in fact, quite 
different types of networks exhibit qualitatively very sim- 
ilar phenomena. That circumstance likely points toward 
the possibility of generalized classification scheme for sta- 
tistical models in a networked world. We will discuss the 
possibilities of such a classification elsewhere Q. 

One such ubiquitous phenomenon is the so-called in- 
verted Berezinskii-Kosterlitz-Thouless transition (BKT) 
0, named after the distinctly weak infinite-order BKT 
phase transition observed for the 0(2)-vector model on a 
two-dimensional lattice 4 6] . In the original BKT transi- 
tion, spin-waves curl into vortex structures, whose pair- 
creation and unbinding dominates the low-temperature 
behavior of the system. Renormalization group analysis 
demonstrates the existence of a line of stable fixed points 
in the low-T regime that is parameterized by the tem- 
perature T, resulting in an equally T-dependent scaling 
behavior and an exponential singularity for the corre- 
lation length there. Analogously, for several Ising spin 
systems on certain networks, a similar T-dcpcndcnt line 
of absorbing fixed points has been found, but inverted 
in the sense that it occurs only in the /wg/i-temperature 
regime, for all T > T c 04H1- Clearly, such a behavior can 
not be related to spin-waves due to the discrete, O(l)- 



symmetry of Ising spins, nor to any topological defects 
such as vortices due to the long-range nature of bonds in 
any of these networks. 

Here, we use an exact real-space renormalization group 
(RG) to study the ferromagnetic Ising model on the re- 
cently introduced set of Hanoi networks [l2| - tl3 | . Hanoi 
networks have been used to study phenomena as diverse 
as diffusion [rJ 0] , synchronization [l2j , the exclusion 
process [HI, percolation quantum transport [l7| . 
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and the vertex cover problem [18(. These hierarchically 
constructed networks possess a regular degree or have 
an exponential degree distribution, akin to the Watts- 
Strogatz Small Worlds [19j . Unlike scale- free networks 
often used to model social phenomena [2C ,|21| , they have 
a more "physically" desirable geometry 22| of a lattice 
backbone mixed with a small-world hierarchy. With re- 
gard to spin models, certain Hanoi networks might com- 
bine finite transition temperatures with susceptibilities 
that exist in the high-temperature regime. They hold the 
potential for an interesting, analytically tractable inter- 
polation from that small-world behavior toward that of 
a finite dimensional lattice 0, [24| by re- weighting long- 
range bonds. As some of those features are lacking in 
other hierarchical networks [1, [25l - l27| , this will serve as 
focus of future investigations. For these Hanoi networks, 
we determine the location of stable and unstable fixed 
points in their phase diagrams and find a range of inter- 
esting RG-flows. Our results demonstrate that a scale- 
free network or highly heterogeneous coupling strengths 
is not a prerequisite for the BKT transition. Moreover, 
we provide analytically tractable one-parameter interpo- 
lations between a ferromagnetic transition governed by 
unstable fixed points and BKT-likc transitions. In gen- 
eral, we find that the critical divergence of the correlation 
length with temperature becomes increasingly singular 
- from power-law over BKT-like to a full essential sin- 
gularity - for increasing relative strength of long-range 
couplings in the network. 

In the following Section, we describe the Hanoi net- 
works. The analysis of the phase diagrams and the RG- 
flow for the Ising ferromagnet on these networks is dis- 
cussed in Sec. Mil In Sec. IIV1 we introduce families of 
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FIG. 1: (Color Online) Depiction of the 3-regular network 
HN3 on a semi-infinite line. 



interpolating networks to reveal a more comprehensive 
set of regimes, each with its own characteristic type of 
phase transition, and we conclude with a discussion of 
our results in their implications in Sec. [V] 



II. GEOMETRY OF THE HANOI NETWORKS 

Each of the Hanoi networks possesses a simple geomet- 
ric backbone, a one-dimensional line of sites n, < n < 
TV = 2 k (k — > oo). Each site is connected to its near- 
est neighbor, ensuring the existence of the ld-backbonc. 
To generate the small- world hierarchy in these networks, 
consider parameterizing any integer n (except for zero) 
uniquely in terms of two other integers (i, j), i > 0, via 



= 2*(2j + l). 



(1) 



Here, i denotes the level in the hierarchy whereas j > 
labels consecutive sites within each hierarchy. For in- 
stance, i = refers to all odd integers, i = 1 to all 
integers once divisible by 2 (i.e., 2, 6, 10,...), and so on. 
Depending on its level of the hierarchy, any site has also 
small-world (i.e. long-range) bonds to more-distant sites 
along the backbone, according to some non-random rule. 
For example, we obtain a 3-regular network HN3 by con- 
necting also 1 to 3, 5 to 7, 9 to 11, etc, for i = 0, next 2 
to 6, 10 to 14, etc, for i = 1, and 4 to 12, 20 to 28, etc., 
for i = 2, and so on, as depicted in Fig. [TJ 

Previously [12], it was found that the average chemical 
path between sites on HN3 scales as 



,HN3 



(2) 



with the distance I along the backbone. In some ways, 
this property is reminiscent of a square-lattice consist- 
ing of TV lattice sites. The diameter (=diagonal) of this 
square is also ~ y/~N. 

While HN3 (and HN4 Q3) are of a fixed, finite de- 
gree, we introduced here convenient generalizations of 
HN3 that lead to new, revealing insights into small- world 
phenomena. First, we can extend HN3 in the following 
manner to obtain a new planar network of average degree 




FIG. 2: (Color Online)Depiction of the planar network HN5, 
comprised of HN3 (black lines) with the addition of further 
long-range bonds (green-shaded lines). Note that sites on the 
lowest level of the hierarchy have degree 3, then degree 5, 7, 
etc, comprising a fraction of 1/2, 1/4, 1/8, etc., of all sites, 
which makes for an average degree 5 in this network. 



5, hence called HN5: In addition to the bonds in HN3, 
in HN5 we also connect all even sites to both nearest 
sites within the same level of the hierarchy i(> 1). The 
resulting network remains planar but now sites have a 
hierarchy-dependent degree, as shown in Fig. [2J To ob- 
tain the average degree, we observe that 1/2 of all sites 
have degree 3, 1/4 has degree 5, 1/8 has degree 7, and 
so on, leading to an exponentially falling degree distribu- 
tion of V {a = 2i + 3} oc 2~\ Then, the total number of 
bonds L in the system of size N = 2 k is 

fc-2 

2L = 2 (2k- l) + ^2(2i + 3)2 k ~ 1 - 1 



i=0 



5x2* 



and thus, the average degree is 

2L 

(a) = — 

V 1 TV 



(3) 



(4) 



In HN5, the end-to-end distance is trivially 1, see 
Fig. [21 Therefore, we define as the diameter the largest of 
the shortest paths possible between any two sites, which 
are typically odd-index sites furthest away from long- 
distance bonds. For the TV = 32 site network depicted in 
Fig. [21 for instance, that diameter is 5, measured between 
site 3 and 19 (starting with n = as the left-most site), 
although there are many other such pairs. It is easy to 
show recursively that this diameter grows as 



,HN5 



= 2 



+ 1 - log 2 TV. 



(5) 



We have checked numerically that the average shortest 
path between any two sites appears to increase logarith- 
mically with system size TV as well. 

The networks HN3 and HN5 have the convenient but 
(from a mean-field perspective) unrealistic restriction of 
being planar. In fact, with a minor extension of the def- 
inition, it is easy to also design Hanoi networks that are 
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FIG. 3: (Color Online)Depiction of the non-planar Hanoi net- 
work HNNP. Again, starting from a ld-backbone (black lines), 
a set of long-range bonds (blue-shaded lines) is added that 
break planarity but maintain the hierarchical pattern set out 
in Eq. (JTJ). The RG on this network remains exact. Note that 
sites on the lowest two levels of the hierarchy have degree 3, 
then degree 5, 7, etc, comprising a fraction of 1/2, 1/4, 1/8, 
etc., of all sites, which makes for an average degree of 4 in 
this network. 



both, non-planar and fully renormalizable. The simplest 
such network, which we dub HNNP, is depicted in Fig. [3] 
To obtain the average degree, we observe that 1/2+1/4 
of all sites have degree 3, 1/8 has degree 5, 1/16 has 
degree 7, and so on, leading to an exponentially falling 
degree distribution, as for HN5. The total number of 
bonds L in the system of size N = 2 k is 

fc-i 

2L = 3 x 2 k - 1 + ^(2i- 1) 2 k ~ l + 2k + (2k - 2), 



= 4 x 2 k - 4, 
and thus, the average degree is 



2L 



4. 



(6) 



(7) 



Here, too, it is easy to see that the shortest paths between 
sites increases logarithmically with system size N. Note 
that in the same manner we extended HN3 to HN5 by 
connecting even-indexed sites within each hierarchy, we 
will discuss such an extension of HNNP toward a non- 
planar network of average degree 6, HN6, in Sec. II V Bl 
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FIG. 4: (Color Online)Depiction of the (exact) RG step for 
the Ising model on HN3. The step consists of tracing out 
odd-labeled variables x n ±i in the top plot and expressing the 
renormalized couplings (L' , K' ) on the bottom in terms of the 
old couplings (Lo, Ko, Ki). Note that the original network in 
Fig- [U does not contain couplings of type (Lo,L' ), but that 
they certainly become relevant during the process. 



A. Ising Ferromagnet on HN3 

The RG consists of recursively tracing out spins level- 
by-level in the hierarchy [12J . In terms of Eq. ([I} , we start 
by tracing out all sites with n odd, i.e. i = 0, then those 
n which are divisible by 2 only once, i.e. i = 1, and so 
on. We can always relabel all sites n after any RG step 
by n — > n/2, so that we trace out the respective odd- 
relabeled sites at any level. It is apparent, for instance 
from Fig. [TJ that odd-labeled sites are connected to their 
even-labeled nearest neighbors on the backbone, say, by 
a coupling K (= (3J ). At any level, each odd-labeled 
site x n ±\ is also connected to one other such site x n ^\ 
across an even-labeled site x n with n = 2 (2j + 1) that 
is exactly once divisible by 2. Let us call that coupling 
K\ (= (3Ji). The basic RG step is depicted in Fig. 0] and 
consists of tracing out the two sites x n ±i neighboring the 
site x n for all j with n = 2 (2j + 1). 

We can section the Ising Hamiltonian 



III. RG FOR THE ISING MODEL 

In this section, we study Ising spin models on HN3, 
HN5 and HNNP with the renormalization group (RG) 
Q. First, we consider RG for the Ising model on HN3. 
This allows us to introduce our procedure although the 
result turns out to be trivial in the sense that there is no 
finite-temperature transition. Especially, it is actually 
identical (except for one bond-type) to the exact RG for 
HN5, and it is almost identical to the treatment below 
for HNNP. 



-m = Yl (~l3H n )+K{K 2 ,K 3 ,...), (8) 

n=l 

where 1Z contains all coupling terms of higher level in the 
hierarchy, and each sectional Hamiltonian is given by 

-/3H n = 41 + L (x n - 2 x„ + x n x n+2 ) + 

K {x n -2X n -l + X n ^lX n + X n X n+ l + X n+ lX n+2 ) 

+K 1 x n - 1 x n+ x, (9) 

where (Kq, Ki, Lq) are the unrenormalized couplings de- 
fined in Fig. 2] and I is a constant that fixes the overall 
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energy scale per spin. (There are effectively 4 spins in- 
volved in each graph- let, as those at each boundary are 
equally shared with neighboring graph-lets.) While cou- 
plings of the type Lq between next-nearest even-labeled 
neighbors emerges that are not part of the network ini- 
tially in HN3, they do emerge during the RG step (other- 
wise the system of recursion equations would not close), 
see Fig. 0] 

To simplify the analysis, we introduce new variables 
similar to inverse "activities" fal, 



C = e" 4/ , K = e~ iK °, A = er iL ", ,i = e" 2 ^ 



(10) 



which ensure that the RG flow only contains algebraic 
functions and, for the ferromagnetic model, remains con- 
fined within the physical domain < K,X,fx < 1. Thus, 
we rewrite Eq. ([9]) as 



-PH n 



c- 



■ ^(x rl -2X n -l+X rl - 1 X n +X n X rl+1 +X rl+1 X rl + 2) 
-lX n+1 



^-i(x n -2X n +X n X n + 2 ) ^—^X n -lX n +l QJ\ 

Tracing out the odd-labeled spins, we have to evaluate 



cr 



E E e- 

i-i=±l} {a;„ + i=±l} 



5 \-\( X -n-2X n +X n X n + 2 ) 



(12) 



K -i(x„_2+2a:„+a:„ + 2) i 



,2) 



+ fJ/K i(Xr l -2-X n + 2) _|_ K ^(x n -2+2x rl +X rl + 2) 
= (C')~^ (A') ^I X r l -2X rl + 2 r K '\ -\(x n -2X n +X n X n + 2 ) 

for the remaining spins in terms of the renormalized 
quantities C , k' , A' . Of the eight possible relations re- 
sulting from the combinations x n -2, x n , x n+ 2 = ±1, only 
three are independent. After some algebra, we extract 
from those the RG recursions: 



k\- 



A' = 
C = 



2(1 + //) 

1 + 2/AK + K 2 ' 

(1 + k) 2 (1 + M ) 
2(1 + 2^k + k 2 ) 

C 2 



(13) 



\/2(l + K)(l + /i) 2 s/l + + K 2 



Note that for couplings in higher levels of the hierarchy 
it is K'i = Ki + i for i > 1; correspondingly, these cou- 
plings, and hence, /z, will not renormalizc. Instead, they 
retain their "bare" value /x 2 determined by the tempera- 
ture, kT/J = — 2/ln/z. In this sense, we will use /z as a 
measure of temperature throughout. 

Only half of the contribution to the renormalized en- 
ergy scale is originating with the sectional Hamiltonian 
in Eq. @, since at the next level two such sections are 
combined into one, making C' cx C 2 . While we do not 
consider the recursions for C in this paper, they are es- 
sential to reconstruct the free energy for each system, 
and will be analyzed elsewhere Q. 




FIG. 5: (Color Online)Phase diagram for the Ising ferromag- 
net on HN3. The stable high-temperature fixed point at 
K* = A* = 1 is marked by a black dot, the line of unsta- 
ble fixed points at low temperature by a thick red line for 
k* = and i < A* < 1. For homogeneous initial conditions 
in Eq. I|14p . the ensuing RG flow for each fi (blue lines) starts 
on the line A^ ' = 1 and evolves to a higher value of re' 1 ' but 
smaller value of A' 1 ' before it rapidly veers off toward the 
high-temperature fixed point. RG flow lines are drawn here 
for equal increments in /i = V «;( > . Note that lines are al- 
lowed to cross, as each flow line is parameterized by a different 
value of /i. 



Eqs. (|13|) provide recursions order-by-order in the RG 
for the evolution of the effective couplings characterizing 
increasingly larger scales of the network. To facilitate this 
RG flow, we need to specify initial conditions for a par- 
ticular physical situation realized in the unrenormalized, 
bare network. Here, we restrict ourselves to networks 
with uniform bonds (although many interestin g ch oices 
are conceivable, such as distance-dependence [81 l23l. |24| ) . 
For HN3 this implies that we chose J = 1 as our energy 
scale, such that Ki = j3J = (3 and I = Lq = initially, 
or in terms of Eq. (|10l) . 



C (0) = A (o) = ^ 



K<°) = M 2 = e" 4/3 . 



(14) 



Kq and V 



Searching for fixed points K = Kq 
L = L* , i.e. k' = n = n* and A' = A = A* in Eqs. (JTc 
immediately provides the trivial, high-temperature solu- 
tion k* = A* = 1, i.e. Kq* = L* = 0. Further analy- 
sis yields only a line of (unstable) strong-coupling fixed 
points, 



0, A* 



1 



/' 



(15) 



extending from A* = | for low temperatures, /i = 0, to 
A* = 1 for T — > oo, where fj, = 1, see Fig. [5] Even at 
T = 0, only the renormalized backbone bonds Kq provide 
strong coupling, the emerging long-range bonds Lq only 
exert limited coupling strength. 

Local analysis near the fixed points with the Ansatz 

K n ~ K*+e„, A„ — A*+<5„, e„,<5 n <l (16) 
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FIG. 6: (Color Online)Depiction of the (exact) RG step for 
the Ising model on HN5. This step is identical to that for HN3 
in Fig. |4]aside from the extra link L\ spanning between x n -2 
and Xn+2 (top), which contributes to the renormalization of 
L' (bottom). 



reveals that the high-temperature fixed point is always 
stable and corrections decay exponentially, where the ex- 
ponential contains a factor of 2? = y/W. At the low- 
temperature line of fixed points in Eq. (|15[) we find 



co(l + A*) 



1 



(17) 



which is divergent for all T > 0, i.e. < /j, < 1, making 
the fixed point at T = unstable. For any fixed point, 
there is no linear expansion possible that would yield 
critical exponents. For the initial conditions in Eqs. (|14[) 
corresponding to uniform couplings throughout the un- 
rcnormalizcd network, the RG flow always evolves to the 
high-temperature fixed point, as Fig. [5] shows. Thus, the 
ferromagnet on this network behaves similar to a Id Ising 
model. 



B. Ising Ferromagnet on HN5 

As shown in Sec. [Til HN5 is basically an extension 
of HN3, created by adding a new layer of links to each 
level of the hierarchy. As is apparent from the foregoing 
discussion in Sec. IIII A[ these additions correspond pre- 
cisely to new renormalizable operators (here, the bonds 
L) that inevitably emerge during the RG of HN3, see 
Fig. |U In HN5, these new operators are simply deemed 
an original feature of the network, hence, maintaining the 
RG as an exact procedure. Consequently, the RG itself 
hardly changes, see Fig. [6j it merely differs by one extra 
link in the graph-let, L\, compared to that for HN3 in 
Fig. 01 In Eq. ©, it only adds the term Lix n ^ 2 x n +2 
to the sectional Hamiltonian and, like Lq itself, L\ does 



not get traced over in the calculation in Eq. (|T2l) . Wc 
can introduce these new bonds as yet another free, non- 
renormalizing coupling in the RG and choose, to wit, 



Li = yKt, i.e. 



-4L, 



2y 

= t L ■ 



(18) 



This merely contributes a factor of yT 4 x ^~2x n+2 ^ ^he 
unprimed side of Eq. (|12p , which correspondingly alters 
only the recursion for A' in Eq. (|T3|) by a factor of /i 2y . 
Otherwise using the same definitions as in Sec. IIII A| wc 
obtain the RG recursions for the Ising ferromagnet on 
HN5:[3i 



A' 

c 



tA 



- ,, 2 v 



2(1 + ^) 

1 + 2flK + K 2 ' 
(1+ K ) 2 (l + At ) 

2(1 + 2^ + k 2 ) 



(19) 



c- 



Kfl 



\/2(l + «)(! + A*) 2 V 1 + 2 M« + 



Accordingly, due to the bare existence of the L\ bond, we 
will have to change the initial conditions from Eq. (fH)) 
to 



C<°> = 1. 



A (0) = fj, 2y . (20) 



Analyzing these recursions for fixed points, n' = K = 
k* and A' = A = A*, wc find that the addition of the ex- 
tra long-range bond has eliminated the high-temperature 
fixed point found in HN3. At low temperatures, we find 
similar to Eq. (fT5)) in HN3 a line of fixed points 



0, A* 



2;/ 



+ 



(21) 



which here extends over the entire domain for the long- 
range bonds, < A* < 1 for < /i < 1. Note that al- 
though y represents a continuous interpolation between 
HN3 and HN5, there is a singular limit at y —¥ toward 
an isolated point corresponding to HN3, see Eqs. (|T5|) . In 
the following, we only treat the case of couplings that are 
homogeneous throughout the unrenormalized network, 
y = 1. Consideration of the rich set of transitions oc- 
curring for the family of networks parameterized by in- 
terpolating < y < 1 is deferred to Sec. IIVI 

Dividing out the k* = 0-solution, further analysis of 
Eqs. for y = 1 reveals yet another line of fixed 

points given by 



1 r 
2 



(1 - fx) n + y/(l + Jij (fi 3 - 3 M 2 + 8{i- 4)J , 

(22) 



A'" 



2 - /i + /i 2 + vOL+70 G" 3 - 3m 2 + 8/i - 4) 



which can be expressed most simply in closed form as 
1 



A* 



k* - 1 + ^5 + 2k* +5(k*) 2 +4(k*) 3 (23) 
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FIG. 7: (Color Online)Phase diagram for the Ising ferro- 
magnet on HN5 with homogeneous bonds, y = 1. Unlike 
for HN3 in Fig. \S\ there is no high-temperature fixed point 
here, instead a line of stable fixed points marked by a thick 
black line that extends from the strong-coupling regime at 
k* = and < A* < ^ = 0.309017 ... all the way to the 
high-temperature regime. For the homogeneous-coupling ini- 
tial conditions in Eq. (|20p . the ensuing RG flow (blue-shaded 
lines) starts on the diagonal (dashed line) and always evolves 
toward lower values of k* and A*. In the low-temperature 
regime, the flow reaches strong coupling, k* — 0, for values 
below [i c = (j)- 1 = 0.618033... (cross at k (0) = A (0) = ill). 
Above that, the flow terminates on the line of finite-coupling 
fixed points, all the way to infinite temperature, /Lt — > 1. 



by eliminating \i. As we will see, these relations lead 
to physical fixed points only within a limited range of 
the temperature [i. We have plotted the phase diagram 
for HN5 at y = 1 in Fig. [71 It also shows the RG flow 
for the initial conditions in Eqs. ([2TJ1) . which starts on 
the diagonal, representing all-equal bonds for the homo- 
geneous network. For these initial conditions, the flow 
always evolves toward smaller values of k and A, i.e. 
stronger coupling. But there is a notable transition where 
the attained fixed point jumps from the low-temperature 
branch in Eq. (|2"T]) characterized by k* = 0, i.e. a solidly 
frozen backbone, to the branch given Eq. on which 
both k* and A* are finite. We obtain this transition 
point by evaluating Eqs. (|22II23[) for k* — 0, which yields 
and a critical temperature of 



A * = ^ = 0.309017 



1 kT c 
4> ° r — 



In 



= 4.156173842. 



/'<■ 



(24) 



where <j> = (^5 + l) /2 = 1.6180339887 ... is the "golden 
ratio" [l3j |. For bare couplings at this temperature, 
marked by a cross in Fig. [71 the RG flow marginally 
reaches the strong-coupling limit. In this network, for 
these initial conditions the RG flow never reaches an un- 
stable fixed point such as the unstable portion of Eq. ([21]) , 
marked by a red-shaded line in Fig. [7J As we will sec be- 
low, this circumstance will change when we weaken the 
impact of long-range couplings. 



0.8 



0.6 



0.4- 



0.2 







0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 

FIG. 8: (Color Online)Plot of the two eigenvalues a± (blue- 
shaded lines) resulting from Eq. (|26|) . In the physical regime, 
4 < H < 1, marked left by the vertical dashed line, both 
eigenvalues satisfy 1 > a+ > > Q- > — 1. Just below 
that regime, a+ exceeds unity, re-enters, and both eigenvalues 
soon (at fj, = 0.6117 . . .) separately disappear into the complex 
plane. 



1. Fixed-Point Stability 

We determine the stability of the fixed points with a 
local analysis using Eq. (|16[) . At y = 1, inserted into the 
strong-coupling solution, Eq. (|21[) , provides 



tn+l 



M 2 (i+^r 



(25) 



and 5 n oc e„. Thus, this line of fixed points is stable 
only for /i smaller than the solution of 1 = /U (1 + //), i.e. 
< /i < 4>~ l ; for larger values of fj, the line is unsta- 
ble. In terms of A*, this implies that the stability of the 
line of strong coupling fixed points changes exactly at 
the location in the phase diagram, Fig. where the line 
of high-temperature fixed points in Eq. (|23[) intersects, 



at ( ^t,0J. Applying the Ansatz on the line of high- 



temperature fixed points in Eqs. 
algebra 



yields after some 



e «+i = -r- [2-M 



1 + /i 



2 - 3/z - n 2 + y/(l + n) (/i 3 - 3/i 2 + 8^ £4) 
2(1-M)M 



/i 1 / [i 3 - 3fi 2 + 8fj, - 4 
2 ~ 2 V 1+Ji ' 



(26) 



The two eigenvalues a± characterizing this system are 
easy to find but lead to messy algebraic expressions; in- 
stead of presenting them, we merely plot their values as 
a function of fi in Fig. [51 It shows that instability and 
a delicate root-singularity (at jj, = 1 + | \/ 6VH4 — 27 — 

5/\/6v / Tl4- 27 = 0.6117 . . .) lurk just outside the phys- 
ical regime, 4>~ 1 < [i < 1, yet inside of that regime the 
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eigenvalues both satisfy \a±\ 
fixed points remains stable. 



< f and the entire line of 



2. Divergence of the Correlation Length at T c 

The effect of the critical point at ll c = \ on the scaling 
behavior of the Ising model deserves special attention. 
For instance, we can associate a correlation length £ at 
a given temperature ll to the system by considering the 
approach to the stable fixed point, which is derived from 
the dominant eigenvalue a+ (fi) < I of Eq. (|26|) affecting 
the asymptotic behavior of the backbone couplings ft„, 

e n ~ ( a +) n e o ~ e"^, 

defining a cut-off scale n* . Aside from a few initial tran- 
sient rescalings, associated with local structure of the 
network, the appearance of the physical state of the sys- 
tem remains scale invariant, e„+i ~ e n , for system sizes 
N = 2™ with n < n* . In turn, for larger systems with 
n > n* , correlated domains of size 



(27) 



become mutually decorrelated to attest to the off-critical 
macroscopic state of the system. Only for /j, — > fx c does 
the correlation length diverge: Expanding the eigenvalue 
resulting from Eqs. (|26[) for ll — > jij , we obtain a + ~ 
1 — a (/x — jU c ) with a = 2\/5. In fact, the correspond- 
ing expansion around the strong-coupling fixed point in 



Eq. (HJ) for fx 



we obtain «j 



1 - a (ji c - fx) 



with the same constant a, so that we can write a+ 
1 — a \ fi — fx c \ and get 



(28) 



This exponential divergence neither resembles the power- 
law divergence for a finite-order phase transition of Ising 
systems on a lattice nor the infinite-order transition char- 
acteristic of BKT that we discuss below. 



C. Ising Ferromagnet on HNNP 

HNNP represents a drastic change in the geometric 
properties over the other Hanoi networks, HN3 and HN5, 
while nonetheless preserving the exact RG. Here, a cross- 
ing set of bonds K\ in the elementary graph-let shown in 
Fig. [5] render this network non-planar. Yet, its average 
degree of 4 is between that for HN3 and HN5, and we will 
see that the effective strength of the renormalized cou- 
pling is intermediate between HN3 and HN5, too: Unlike 
HN3, it does have an ordered low-temperature regime, 
but unlike for the transition in HN5, its high-temperature 
regime above T c only possesses a high-temperature fixed 
point of vanishing coupling strength; not enough to sus- 
tain the partially ordered state found in HN5. 



K, K] 



/ K„ / K„ K„ \ K \ 
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n-2 n n+2 



FIG. 9: (Color Online)Depiction of RG step for the Ising 
model on HNNP. This step differs from that for HN3 in Fig. H 
and HN5 in Fig. [6] as there are a set of crossing long-range 
bonds in each elementary graph- let before the RG step (left). 
But the resulting graph-let after the RG step is again identical 
to that for HN3 and HN5. 



The structural feature of the crossing bonds is reflected 
in the Ising Hamiltonian, of course, and we find for its 
sectional Hamiltonian 

-/3Hn = 4/ + L Q (x n - 2 x„ + x n x n+2 ) + 
K (x n - 2 x n -i + x n -ix n + x n x n+ i + x n+ ix n+2 ) 
+Ki (x n _ 2 x n +i + x„„ia;„ + 2) (29) 



Following the same procedure as in Eqs. (jlOllf 3[) yields 
the RG recursions 



«; A 



(i + m) 2 

{1 + Llnf 



A' 



{1 + uikY 

c = c 2 



(30) 



(I + /i) (k + h) (I + lik) 



with the same initial conditions as for HN3 in Eq. (JTJ 

Once more inspecting these recursions for fixed points, 
k' = k = k* and A' = A = A* , we find a strong-coupling 
line of fixed points at 



k* =0, A* = ii , 



(31) 



which extends over the entire domain for the long-range 
bonds, 0<A* < lfor0</j< 1. We can similarly 
identify a high-temperature fixed point at 

«* = A* = I (32) 

that is valid for all fi. Further analysis of Eqs. (|3U|) reveals 
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yet another line of fixed points given by 

rf = I^f< A* = ii^£ (33) 

It is clear that this line of fixed points is physical merely 
in the range ^ < /j, < ^ beyond which one or both of k* 
and A* leaves the unit interval. We can express this line 
in closed form, 

k* = A* + \/A* - 1, (34) 

by eliminating /i. 

We have plotted the phase diagram for HNNP in 
Fig. [10] It also shows the RG flow for the initial condi- 
tions in Eqs. (fl"4"|) . which starts on the right-hand vertical 
axis, A^ ' = f . For these initial conditions, there appears 
to be an ordinary, finite-order phase transition, with the 
RG-flow diverging from an unstable fixed point toward a 
solidly frozen ordered regime at lower temperatures and 
a plain disordered regime above. Yet, the transition is 
non-universal, for instance, dependent on the initial ra- 
tio between Kg and A'i bonds. There is no patchy order 
in this system, as there are apparently not enough long- 
range bonds to enforce it, compared to HN5. The actual 
critical temperature T c is the transcendental solution of 
an infinite integration of the RG-recursions in Eq. (f3"0)) 
from the initial conditions in Eq. (|I4[) . which can be ob- 
tained to any accuracy by a shooting procedure. With 
that, we obtain 




0.2 0.4 0.6 0.8 1 



FIG. 10: (Color Online)Phase diagram for the Ising ferro- 
magnet on HNNP. Similar to HN3 in Fig. \5\ there is again a 
high-temperature fixed point here (black dot). Though, more 
akin to HN5 in Fig. a line of fixed points extends along 
k* = and < A* < 1, which for A* < <f>~ 2 is stable (thick 
black line) and unstable above (thick red line). Both regimes 
are connected by a (red) line of unstable fixed points given by 
Eq. (O with an intercept at A = <?T 2 = 0.381966 ... For the 
homogeneous-bond initial conditions in Eq. (|14[) . the ensu- 
ing RG flow (blue lines) starts from A' ' = 1, evolves toward 
the diagonal k"' = A' 1 ' first, and proceeds in an oscillatory 
manner. For initial temperatures below the critical point 
= 0.30671... (x), the RG flow eventually veers toward 
strong coupling, k* = 0, and above, to the high-temperature 
fixed point. The particular fixed point governing the flow for 
these initial conditions is marked by a red dot. 



H c = 0.553814426157623 . 



3.3845207164300986. 



(35) 



(36) 



marked by a blue cross in Fig.QJJ Inserted into Eqs. ([33| . 
this places the unstable fixed point governing the crit- 
ical behavior for these particular initial conditions at 
(k*,A*) = (0.4547454105..., 0.64908641578...), labeled 
by a red dot. 



1. Fixed-Point Stability 

We explore the stability of these fixed points with a 
local analysis, using again the Ansatz in Eq. (fTB"]) . First, 
we consider the strong-coupling fixed point in Eq. (1311) . 
Expanding the RG-recursion for k' with the Ansatz yields 



e„ + i = /j 2 (1 + /i) 2 e„, 



(37) 



immediately proving that this fixed point is stable only 
for /x < -x , where the n = 0-line intercepts with the other 
line of fixed points in Eqs. (|3"3")l or (|34[) . For values above 
that, = becomes unstable. Conversely, expanding 



around the high-temperature fixed point in Eq. 

1-M 



£n+l 
5n+l 



1 + M 
2(1-M) 



l + pt 

which possesses the eigenvalues 



gives 



(38) 



a± 



1-M 
2 (1 + At) 



1± a/1 +8 



1 + At 



l-/x 



(39) 



The lower branch for all < [i < 1 remains confined to 
— 1 < a_ < 0, while 1 < a + < 2 for /i < i; only for /i > i 
is this fixed point stable. In the overlap, \ < [i < 4, of 
the stable regimes of the two previous fixed points, we 
find the remaining line of fixed points in Eqs. (|33p . which 
prove to be unstable. To show that, we insert the Ansatz 
to obtain the system 



Cn+l 
On+l 



fl 2 + 2^1 - 1 1 - fl - fl 2 



I -ll 2 



Sn, 



1+fJ, 



(40) 
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The eigenvalues a± of this system are given by 

/i 2 + 2/i - 1 ± Jl + Afi + 2^ 2 - 12^ 3 - 7^ 

a± = w^W) m 

for which cv+ > 1 and > cv_ > — 1 on the physically 
relevant interval i < /i < i. Thus, any fixed point 
on that line has an unstable direction and a stable but 
oscillatory direction, as is apparent for the flow in Fig. If 01 

2. Divergence of the Correlation Length at T c 

Although we have found that the phase diagram for 
HNNP consists of lines of stable and unstable fixed 
points, only one specific unstable fixed point controls the 
scaling behavior for any specific incarnation of the sys- 
tem. Any such incarnation is selected by the initial choice 
of couplings A 1 - - 1 ) as a function of temperature //. 

The controlling fixed point in Eq. (|33[) is a non-trivial 
consequence of these initial conditions. It is in particu- 
lar a function of the specific critical temperature \x c for 
which it is reached. The induced scaling behavior for 
nearby temperatures is thus non-universal. Yet, for a 
given /x c , the local analysis for such an unstable fixed 
point proceeds identically to that for Ising systems on a 
lattice @. Here we have two scaling fields, whose fixed- 
point analysis leads to two eigenvalues in Eq. ([4"TT) that 
yields one relevant and one irrelevant variable. From the 
eigenvalue a + (> I) for the relevant variable, we obtain 
for the correlation-length exponent 

v= — = ; — v « 10.8301 .... (42) 
V+ lna+(ji c ) 

such that 

£~|/i- Mc p. (43) 

Thus, we find that the correlation length on approach to 
the fixed point merely diverges with a power law here. In 
this particular case, the exponent v happens to be very 
large, such that it would be difficult to distinguish this 
divergence from the exponential kind numerically. 



y=0.2 

K 




v=0.4 



K 




IV. INTERPOLATION BETWEEN HANOI 
NETWORKS 

We have already observed in the construction of HN5 
in Sec. MI Bl that it is easy to promote the L-couplings, 
which inevitably emerge during the RG, to be associated 
with an actual bond in the network. Here, we will fully 
exploit this fact to obtain a one-parameter family of prob- 
lems with various regimes of phase behaviors, based on 
the interpolating parameter y in Eq. (|18[) . This procedure 
has also been used in Ref. 17J and a similar construction 
can already be found in Ref. In particular, we dis- 
cover transitions between such regimes as a function of 



FIG. 11: (Color Online)Plot of k* in Eq. (J44J) as a function of 
fi 2 for y as defined in Eq. (|18p for the interpolation between 
HN3 and HN5. The generic cases are represented by y = 
0.2 (top), y = 0.4 (middle), and the special value of y — 
'"^ 2 2) = 0.58 . . . (bottom) at which the branch point (BP) 
"sunsets" out of the physical regime and the behavior becomes 
similar to that for y — 1 shown in Fig. [7] In all cases, the 
dashed line indicates the initial conditions (IC) for the RG 
flow in Eqs. (|19H20[> . k {0) = /i 2 . For n fixed, the RG flow must 



proceed vertically, either up or down, to the nearest stable 
line of fixed points. For low y (top), the IC cross the unstable 
branch below BP which then can not be reached by the flow. 
Once the IC cross above BP (middle) , the flow must pass BP, 
unless BP sunsets (bottom). 
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the parameter that will allow us to clarify the connec- 
tions between the diverse set of behaviors that we have 
discovered in the previous section. 



A. Interpolation between HN3 and HN5 



for the strong-coupling solution, Eq. (|2Tj) . yields 



In Sec. MI B[ we argued for the introduction of small- 
world bonds with couplings Li and developed the RG re- 
cursions in (|19j) assuming a relative strength of these cou- 
plings to those germane to HN3 of the form in Eq. (fT8"|) . 
Here, we will now consider the behavior that results from 
varying the strength parameter y between the two ex- 
tremes already explored, y = for HN3 in Sec. IIII Al and 
y = 1 for HN5 in Sec. iHTBl 

Inspecting these recursions in Eq. (|19[) for fixed points, 
we already found the low-temperature fixed point line in 
Eq. {23}. Aside from this n* = 0-solution, Eqs. (JTSJ) 
further reveal a line of fixed points given by 



= - 
2 



)J? (1 + //) - 2(1 ± y/Vyfa) 



(44) 



2 (1 - fi) + n y (1 + IJ.) ± yjv y (fi) 



abbreviating the discriminant 

V y ( /U ) = (l + /x) [fi 2v (l + n) -4(1 -fi v )(l -//)]. 

(45) 

For y — > 0, this solution morphs into the high- 
temperature fixed point of HN3, see Fig. [SJ But for any 
finite y, these lines of fixed points are non-trivial func- 
tions of n, as depicted for k* (/i 2 ) in Fig. [TT] The dom- 
inant feature in these plots is the root-singularity in k* 
with a branch point separating the upper stable and lower 
unstable line of fixed points. Essentially, three distinct 
generic regimes can be discerned: (1) If the branch-point 
happens to lie below the physical regime, we observe a 
phase transition without access to any unstable point (see 
bottom of Fig. fTTjk a critical point akin to that for HN5 at 
y = 1 in Sec. IIII Bl arises. If the branch point rises into the 
physical regime, here for y < y c = Ig^ZS — 0.584963 . . ., 
then depending on whether the initial conditions of the 
RG flow cross the critical line below or above the branch 
point, we find (2) a transition seemingly of finite-order on 
intercepting the unstable lower branch (see top of Fig. fTTj) 
for which the RG flow never accesses the branch point 
singularity. If, in turn, the initial conditions cross above, 
(3) a BKT-likc transition results because the RG flow 
now must pass the singularity (see middle of Fig. 1111) . as 
we will show below. 



1. Fixed-Point Stability 

As before, we determine the stability of the fixed points 
with a local analysis using Eq. (|16l) . Inserting the Ansatz 



M 2y (l + /^f 



(46) 



and S n oc e„. Thus, this line of fixed points is stable only 
for fi smaller than the solution of 1 = fi v (1 + /i); for 
larger values the line is unstable. In each plot of Fig. [TTl 
that value of /z corresponds to the location on the k = 0- 
line where the non-trivial fixed-point line intercepts. 

Applying the Ansatz on the line of high-temperature 
fixed points in Eqs. (|44|) obtains now 



(1 + aQ - 2) ± 
2(1 + ii) fiv 



+ 2{l-n)nv 

Sn+1 - 2{TTJ) e - (47) 

It is tedious but straightforward to show that the two 
eigenvalues 0^2 characterizing this system indicate an 
upper stable branch for "+" and a lower unstable branch 
for "-" . Yet, the decisive question for the observed phase 
transition concerns the intercept with the initial condi- 
tions of the RG flow, i.e. the relative strengths of bare 
couplings in the network. If the intercept lies below 
the branch-point singularity, the phase transition is sim- 
ilar to that discussed in Sec. IIII C[ although here the 
high-temperature phase flows not into a single point but 
rather into a temperature-dependent line of couplings, 
similar to (but inverted from) BKT. Apparently, such 
temperature-dependent couplings in themselves do not 
imply an infinite-order transition. 

In turn, once those initial conditions pass above the 
branch point, such as in the middle panel of Fig. 1111 
there is no stable fixed point above that and the flow has 
to move toward lower values of K during the RG flow. 
Whether there is any intercept with the upper branch 
makes little difference to the generic case here, as the 
critical behavior is controlled now by the flow in the im- 
mediate neighborhood of the root-singularity where both 
branches pinch off. Only in this case do we obtain an 
infinite-order transition of the BKT type. 



2. Divergence of the Correlation Length at T c 

We have already considered the exponential divergence 
characteristic of HN5 for y > y c for the special case y = 1 
in Sec. IIII B 21 Similarly, we have found in Sec. IIII C 21 
that HNNP at y = 1 has a critical point very much like 
that found for HN5 at small y > 0, aside from a difference 
in the high-temperature coupling. For both cases, shown 
in the bottom and top panel of Fig. 111! respectively, fur- 
ther analysis would not alter the qualitative behavior in 
the divergence of the correlation length £. 

Therefore, we restrict ourselves here to the novel case 
where the RG flow must pass near the branch-point sin- 
gularity for a given choice of initial conditions displayed 
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in the middle panel of Fig. [TT] for HN5 with y just below 
y c . As has been argued for a similar case in a scale- free 
version of the hierarchical lattice [|[, for temperatures n 
approaching the critical temperature ii c ever closer from 
below, the flow spends an ever larger number of itera- 
tions near the singularity. This characteristic number of 
RG iterations n* (li) diverges for /i — > /i~ and can be in- 
terpreted as a correlation measure, just as in Scc. lIIIB~2l 
But unlike the n* ~ l/|/z c — ji\ scaling there, we find 
here that n* ~ 1/ \/\ii c — iAi leading the the BKT-likc 
transition also found in many other networks. 

To derive n* , we follow Ref. Q in ignoring the transient 
number of RG steps that would take us from the initial 
conditions to the singularity. Instead, we imagine that 
our flow would start right at the singularity but with an 
infinitcsimally small shift down in temperature, 



dropping terms cubic in e, S and any smaller terms in r. 
All coefficients are algebraic functions of li c and y. The 
proper scaling field u n for the impending analysis is iden- 
tified by multiplying the (5-equation by C and combining 
it with that of e: 



(A + CF)t + (B + CG) e n + CS n 
+ (D + CH)e 2 n + Ee n S n . 



(53) 



The subtlety of this procedure reveals itself now in the 
fact that B + CG = 1 for any ll c and y, and we can write 
with some new (positive) coefficients 



u' (n) 



-Pt 



T = /i c - LI, 



(48) 



so as to dislodge the flow from what is a fixed point for 
/j = jU c . The system then decorrelates, i.e. the RG flow 
escapes the singularity by a finite amount, for sizes N ~ 
2™, i.e. after n > n* iterations of the RG. 

As the transition is obtained right at the branch point, 
the critical temperature \x c derives from 

V y ifx c ) = (49) 

in Eqs. (|44H45I) . We define the critical-point couplings 



using a continuum approach justified for large n. In the 
limit r — > 0, this equation describes a singular boundary 
layer problem easily analyzed with the techniques out- 
lined in Ref. [28|. Instead, we note that the differential 
equation can be solved exactly to give (with uq =0): 




r tan I y PQrn 



(54) 



Thus, the solution exits the boundary layer near the 
branch point (moving down toward the strong-coupling 
fixed point) and becomes finite for y 'PQrn* — > ?r/2, and 
we identify the scaling 







K* (fl c ) 




-At c + ^ (1 + Mc) 


n* 






_ A* (Mc) _ 




_4(1- Mc )+^(1+ Mc )_ 





1 



1 



/' 



(55) 



(50) 

from Eq. (j44|). Such a phase transition only exists for 
both, k c and A c , in the physical regime and, hence, it 
disappears when k c = is reached. The simultaneous 
solution of k c = and Eq. (|49]) yield 



which by Eq. (|27p leads to the divergence in the correla- 
tion length characteristic of BKT, 



e vmJ 



(56) 



Inj 
In 2 



at 



1 

2' 



(51) 



see Fig. [TT] (bottom), with a transition similar to that at 
y = 1 for all y > y c . 

For y < y c , we expand Eqs. (|19[) around the branch 
point to 2nd order in e„, S n with 



Clearly, the physical origin of this singularity is not even 
remotely related to an actual BKT transition. In fact, 
instead of its rarity, confined to very particular lattice 
models, we may find it to be one of a few generic types 
of transition in networks. 



B. Interpolation between HNNP and HN6 









An 




x c + s n 



and eo = So = 0. Note that neither k c nor A c depend on 
Li; it is precisely the explicit appearance of fj, in Eqs. (fT9|) 
that allows a further lst-ovdev expansion in r according 
to Eq. (|48j) . For general y, the solutions of Eq. (|49|) arc 
too messy to display the results, but the nature of the 
recursions is as follows: 

e„+i - At + Be n + C5 n + De 2 n + Ee n 5 n , 

S n+1 ~ FT + Ge n + He 2 n , (52) 



Similar to the step from HN3 to HN5, we can add the 
emerging bonds Li again as an original feature of the 
network, which we parameterize as in Eq. (|18|). As a 
result, we obtain a network that has an average degree 
of 6 with an exponential distribution of bonds which we 
will call HN6, accordingly. While difficult to display, HN6 
consists of the addition of the green-shaded bonds from 
HN5 in Fig. Uto HNNP in Fig. [3] 

As in Eq. (p~8|) . we add the extra L-bonds with an ini- 
tial coupling strength that is a y-fraction of the original 
HNNP bonds; for y —> we approach the pure HNNP 
from before, and for y = 1 we attain HN6 with all equal 
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couplings. In this calculation, we get similar to Eqs. (fT9]) and ([30 
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FIG. 12: (Color Online)Plot of k* in Eq. (J59j| as a function of 
ji 2 for y as defined in Eq. (|18p for the interpolation between 
HNNP and HN6, similar to Fig. 1111 The generic cases are 
represented by y ~ (top), y = 0.3 (middle), and the special 

1h(i+n/2) _ 



value of y — 2 



= 0.54 . . . (bottom) at which 



the branch point "sunsets" out of the physical regime, see 
discussion in Fig. [TT] 



A' 



kX 



; 2 'i 



(« + /*) 2 

(1 +/iK) 2 



(57) 



C" = C 2 



K(l 



(i + m) (« + a») (i + a»«)' 



which is a simple extension of Eqs. (I3TJ1) . to which these 
RG recursions reduce for y = 0. Similar to the interpola- 
tion between HN3 and HN5, this interpolation between 
HNNP and HN6 also possesses interesting phase behav- 
ior. In particular, we also find a non-trivial value of y c 
separating a HNNP-like phase diagram for low y from 
one similar to HN5 for larger y. 

First, we find from Eqs. (|57[) the strong-coupling fixed- 
point line 



0, 



A?: 



/' 



2+2y 



(58) 



which is stable for all (i below the solution of 1 = 
H 1+y (1 + fj,). The remaining fixed points in the physi- 
cal regime satisfy 



1 

M 



1 



2(i 



2~2y 



l±jA y ((i) 



where we abbreviate the discriminant 



A y ( M ) = i-V- w (i-aO 



(59) 
(60) 

(61) 



Note that for y — > 0, the upper branches condense into 
the high-temperature fixed point of HNNP, k* = X* —> 1, 
and the lower branch becomes the unstable fixed-point 
line of HNNP. It is similarly easy to show that, indeed, 
for all < y < 1, "+" corresponds to stable fixed points 
and "— " to unstable ones. 

At this point, it is obvious that the behavior for this 
system parallels that for HN5 above. In Fig. [T^l we dis- 
play for representative values of y the dependence of k* 
in Eq. (|59p as a function of temperature (l, which in each 
panel is equivalent to those in Fig. [TTJ For the small- 
est values of y (here, including also y = 0), the initial 
conditions in Eq. (fT4)) result in an RG flow that avoids 
the branch-point singularity and instead cross a line of 
unstable fixed points, one of which dominates the criti- 
cal behavior that leads to a non-universal power-law di- 
vergence in the correlation length £ as in Eq. (|43[) for 
y = 0. For stronger long-range couplings in the network, 
beyond some value of y for those initial conditions, that 
intercept moves above the branch point, forcing the flow 
through it at the critical point, which now develops a 
BKT-like essential singularity in £. Increasing such cou- 
pling strength further moves the branch point below the 
physical regime, out of reach from any initial condition. 



13 



Here, this requires both, A y (/x) in Eq. (f6"Tj) and k* 
Eq. ([55)1 , to vanish, which occurs for 



in same manner a BKT-likc divergence of £ as in Eq. (|56[) 



/Lt c 



1 

71' 

In [ylt c (1 + /i c )] 

ln/j c 



(62) 



= 0.5431066. 



Hence, for all y > y c , there is no more high-temperature 
line of unstable fixed points and HN6 behaves identical 
to HN5 in Sec. HHBl 

The analysis of the critical behavior for the regimes of 
Uc < y < 1 and for y just below y c proceeds identically to 
Sec. IHI Bl and Sec. II V Al respectively. For the former, we 
just mention that the expansion along the upper stable 
branch obtains 



£n+l 



At 2 - y/Ay (//) 



H (l - y/A y (/*) 



2(1 -M)' 



^n+l — 



1 



l-SjAy (fl) 



(63) 



For example, at y = 1, the behavior of the eigenvalues a± 
of this system qualitatively resemble that in Fig. |8] The 
larger eigenvalue destabilizes (i.e. exceeds unity) when 
the temperature is lowered to 



He 




1 



0.754878 . 



(64) 

or kT c /J = 7.11239..., below which both eigenval- 
ues disappear separately into the complex plane at 
Near fx c , we find for the dominant eigen- 

~ 1 - o|m-Mc| s 



A* = 
value 



3 
4 • 

a j 



with a constant 



1 + "a/(623 - 75 /2 + a (623 - 75^) /2 



2 
3 

6.43855 . . ., leading to the same conclusion as in 
Sec. IIII B~2l with a correlation length given by Eq. (|28| . 

In turn, for values of y just below y c , we locate the 
branch point via A y (p, c ) = to obtain the critical tem- 
perature Hc[y) and find 







K* {fJL c ) 






_ A* (Mc) _ 



(65) 



where we have eliminate all explicit reference to y. Fol- 
lowing the steps subsequent to Eq. (|50|) . we find in the 



V. CONCLUSION 

We have analyzed the fixed-point structure of an Ising 
ferromagnet on a set of Hanoi networks with an ex- 
act real-space renormalization group. Using interpolat- 
ing families of such networks, with the relative coupling 
strength y between backbone and small-world bonds 
as the interpolation parameter, we reveal a number of 
regimes with distinct critical behaviors. While in each 
such regime the critical transition has non-universal fea- 
tures, the characteristics of the transition in each one 
has generic, robust features. For increasing strength, 
we observe that the divergence in the correlation length 
changes from power-law x~ u , to a BKT-like essential sin- 
gularity e 1 /^, then to a plain essential singularity e 1 / 1 , 
on approach to the critical point x ~ \fi c — n\ — > 0. We 
trace the changes from one regime to the next in terms 
of the analytic structure of the RG flow. Finding an enu- 
merable range of such characteristics hints at a possible 
classification of critical behavior of statistical models in 
networks generally, a task we will explore in more detail 
elsewhere Q • For example, critical properties of the kind 
found here have also been observed in community forma- 
tion H and in percolation [ll], EE l29j - l3lj . The existence 
of entire regimes that exhibit essential singularities in 
the divergence of the correlations, as we have found here, 
might explain the surprising prevalence of typically quite 
rare BKT-like transitions in otherwise unrelated network 
models 0. 

There are a number of possible extensions of this work, 
which we are currently exploring. Here, we have merely 
focused on the fixed-point structure and exclusively con- 
sidered the divergence of the correlation length. We are 
currently developing the RG scheme to investigate more 
complicated observables to extract insights into the phys- 
ical origin of these peculiar critical phenomena. Merely 
the addition of an external field as a generator for mag- 
netic properties already makes the analysis much more 
complicated by breaking the ^-symmetry, thus allowing 
for yet another renormalizablc bond linking three spins 
simultaneously [lH . As well, the study of other models 
spin in this context seems desirable. The absence of a bi- 
partite structure (i.e. the existence of odd-length loops) 
in the Hanoi networks promises interesting effects due to 
frustration, which can be calculated with similar rigor for 
an anti-ferromagnetic system (and, at least, numerically 
exact for a spin glass) . The results of those investigations 
are forthcoming. 
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